Effects of the Dzyaloshinskii-Moriya interaction on low energy magnetic excitations in 

copper benzoate 
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We have investigated the physical effects of the Dzyaloshinskii-Moriya (DM) interaction in copper 
benzoate. In the low field limit, the spin gap is found to vary as H'^'^\u^'^{J/ij.bHs) [Hs-. an effective 
staggered field induced by the external field H) in agreement with the prediction of conformal field 
theory, while the staggered magnetization varies as H^^^ and the In^^^ {J/ ^bHs) correction predicted 
by conformal field theory is not confirmed. The linear scaling behavior between the momentum shift 
and the magnetization is broken. We have determined the couphng constant of the DM interaction 
and have given a complete quantitative account for the field dependence of the spin gaps along 
all three principal axes, without resorting to additional interactions like interchain coupling. A 
crossover to strong applied field behavior is predicted for further experimental verification. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.25.+Z 



Recently some novel magnetic properties were discov- 
ered in a variety of quasi-one dimensional materials, such 
as copper benzoate Cu(C6D5COO)23D20[l 0, Yb4As3 
0, 0^ 1^ and BaCu2Si2O70. In these materials, the 
Dzyaloshinskii-Moriya (DM) interaction[3,@,0| plays an 
important role, especially in an applied magnetic field. 
This has stimulated extensive investigation on the physi- 
cal properties of the DM interaction. However, this inter- 
action is rather difficult to handle analytically, which has 
brought much uncertainty in the interpretation of exper- 
imental data and has limited our understanding of many 
interesting quantum phenomena of low-dimensional mag- 
netic materials. 

For copper benzoate, Dender et a?P,|3| found that the 
spin excitation gap shows a peculiar field dependence, 
A ~ in low fields. On the contrary, excitations 

remain gapless in the S=l/2 Heisenberg model below a 
critical field. Oshikawa and Affleck (OA) suggested that 
this field dependence of the gap is due to a staggered mag- 
netic field induced by the DM interaction in addition to 
the staggered g-factor in a uniform fieldjlE^- 

However, 

a satisfactory explanation for the field-dependence of the 
energy gaps in all three directions is still lacking flllll^ . 
It was argued that the inconsistency between the exper- 
imental data and theoretical results might be due to the 
neglect of the interchain coupling and/or anisotropic in- 
teraction terms in the low-field effective model used by 
Oshikawa and Affleck . We believe this issue can 

be clarified by a thorough study of the DM interaction 
and a direct comparison with experiments. 

Copper benzoate is a quasi-lD spin-1/2 antiferromag- 
netic Heisenberg system. The chain direction is the c- 
axis. It contains two types of alternating and slightly 
tilted CuOs octahedra. This leads to two i ncq uivalcnt 
Cu++ ions and an alternating DM coupling [ijj. In an 
applied field, copper benzoate can be modeled by the 



following Hamiltonian, 

i 



(1) 



where the three terms in the summation are the an- 
tiferromagnetic Heisenberg, DM and Zeeman splitting 
interactions, respectively. The exchange coupling con- 
stant J, determined from the neutron scattering mea- 
surements, is about 1.57meV. The DM interaction is 
much weaker than the Heisenberg term. The D-vector, 
primarily aligned along the a" axis, will be determined 
numerically, g" and are the uniform and staggered 
components of the alternating g-tensor as given in Ref. 
for copper benzoate. 

The DM term can in principle be eliminated by per- 
forming a spin rotation about the D vector by an angle 
a — ±itan~^(D/J) on the alternating sites. This results 
in a small exchange anisotropy and an effective staggered 
magnetic field in addition to the g'^-term in The total 
effective staggered field is approximately given by 



g^ - X g" 



H 



(2) 



up to the leading order in D/J. In Ref. 0, OA studied the 
isotropic Heisenberg model with a staggered field after 
eliminating the DM term and neglecting all anisotropic 
terms. They estimated that the DM coupling constant 
is about D = (0.13, 0.0, 0.02) J from the specific heat, 
neutron scattering and ESR measurement dataP,0. 

In order to explore the low-energy properties of copper 
benzoate, we use density matrix renormalization group 
(DMRG) Hi, El to study dhectly the Hamiltonian de- 
fined by Eq. (1). Since neither the total spin nor its 
z-component is a good quantum number, the calculation 



2 



is computer time consuming. Open boundary conditions 
are used and up to 400 states are retained in our calcu- 
lation. The convergence of the results is systematically 
examined and the truncation error is less than 10~^. To 
obtain the values of gaps and magnetizations in the ther- 
modynamic limit, an extrapolation from the DMRG cal- 
culations up to 1000 sites is done. 
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FIG. 1: Spin gaps versus 1/L in different applied fields for 
D = (0.13, 0.0, 0.02) J. The exact diagonalization results (o) 
for small lattices are also shown for comparison. 

Figure ^ shows the spin excitation gap as a function 
of inverse lattice size L in different applied fields H . In 
zero field, there is no gap in the excitation spectrum in 
the thermodynamic limit and the gap decreases mono- 
tonically with L increasing. However, in finite fields, the 
size dependence of the gap becomes more complicated 
and changes dramatically with H increasing. The oscil- 
lations of the gaps are due to the competition between 
the uniform and staggered magnetic fields. The most pro- 
nounced oscillation occurs when the contributions to the 
Zeeman energy from these two fields become compara- 
ble, and the characteristic length scale of the oscillation 
is inversely proportional to the energy gap. Since the 
gap along a"-axis is much smaller than other directions, 
the oscillation looks stronger for iJ ~ (7, 0,0)r than for 
other cases. 

Figure |21 shows the field dependence of the gap with 
D = (0.13, 0.0, 0.02) J. In low fields, we find that the gap 
is well described by the following equation 



(3) 



as predicted by Oshikawa and Affleck. By fitting the 
numerical data with the above equation, wc find that 
the coefficient C" is given by C = (0.0072, 0.059, 0.097) 
for D = (0.13, 0.0, 0.02) J. In this system, the gap is fi- 
nite even in very low field. This is completely different 
from the pure Heisenberg model whose excitations re- 
main gapless until the spin polarization induced by the 
applied field becomes saturated. This difference results 
from the effective staggered magnetic field induced by 




FIG. 2: Field dependence of the spin excitation gap. D = 
(0.13, 0.0, 0.02) J. The solid lines are the fitting curves of the 
numerical results with Eq. the low fields. 



the alternating g-tensor and the DM interaction. Since 
this staggered field couples directly to the strongest spin 
fluctuation mode at g = tt, it gives rise to the instability 
of the ground state and opens a gap in the excitation 
spectrum. In high fields, the gap varies almost linearly 
with H and its slope is proportional to the correspond- 
ing value In the intermediate field regime, the gap 
varies non-monotonically with H and shows a local min- 
imum around H ^ 25T in all three directions. This 
non-monotonic behavior of the gap is also due to the 
competition between the applied field and the induced 
staggered field. However, we find that this competition 
does not lead to a phase transition in the ground state 
in this crossover regime, since the ground state energy 
drops continuously and smoothly with H. 
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FIG. 3: Comparison between numerical results (o and o) and 
the specific heat Cv {+, <) and ESR (>) data for the field 
dependence of energy gaps. Solid lines are the fitting curves 
of the DMRG results to = CgMRG-ff^^^ ln^'''^(J/MB-ffs) 
for D = (0.11, 0.0, 0.0) J with Cdmrg = (0.025,0.048,0.089). 

Now let us turn to the issue whether the model ^ 
is sufficient to describe the low-energy magnetic prop- 
erties of copper benzoate. In the effective field the- 
ory, the critical exponents of the gaps are determined 
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by the relevant and marginal operators in the low field 
limit. However, as this theory contains only the con- 
tribution of the dominant effective staggered field and 
ignores all other anisotropic terms after eliminating the 
DM term, it cannot explain the field-dependence of the 
gaps along all three principal axes. Moreover, since 
the energy gap is very sensitive to the DM interac- 
tion, the value of D-vector determined from this theory, 
D — (0.13, 0.0, 0.02) J0, would also be inaccurate. 

To determine accurately the value of D, we have 
evaluated and analyzed the field dependence of the 
gaps around D = (0.13, 0.0, 0.02) J. By comparison 
with experiments 0, we find that the gaps for D — 
(0.11, 0.0, 0.0) J give the best fit to the experimental data 
in all three directions (Figure 3). It also agrees with the 
recent ESR data for H//c"^' as well as the results of 
Essler for H//b and in low fields for H//c"^. For this 
D value, our numerical results for the gap ratios A^" : 
Ab : Ac" in low fields are 1 : 1.92 : 3.56, in agreement 
with experiments. However, for D = (0.13, 0.0, 0.02) J, 
the low field gap ratios are 1 : 8.2 : 13.5 and the theoret- 
ical results agree with the experimental data only along 
the c"-axis (Figure Therefore the Hamiltonian 
is indeed an appropriate model for copper benzoate and 
the discrepancy between the effective field theory and ex- 
periments is due to the neglect of the anisotropic terms 
rather than due to the interchain coupling as previously 
suggested [sj. 

For the S=l/2 Heisenberg model, the momentum of 
the antiferromagnetic soft mode is shifted from tt to TT±Sq 
in an applied field. This leads to an incommensurate peak 
in the longitudinal structure factor. The momentum shift 
is predicted to be proportional to the magnetization: 

Sq = 2ttM{H). (4) 

This linear relationship between the momentum shift and 
the magnetization was examined by Dender at al for cop- 
per benzoateQ- They found that the momentum shifts 
at H = 3.5, 5, 7T for H//fo were consistent with the the- 
oretical results for the pure Heisenberg model with 
g = 2.059 and J = 1.57meV. However, as the excitation 
spectrum of copper benzoate is fully gaped even in an ar- 
bitrarily small field, unlike in the pure Heisenberg model, 
this issue needs also to be reexamined. To do this, we 
have calculated the spin-spin correlation functions 

The momentum shift Sq can be determined from the peak 
position of S^". From our calculations, we find that 5^" 
behaves differently along the three principal axes. S^^ 
shows a weak incommensurate peak. However the peaks 
of and ° are pinned at g = tt. With H in- 

creasing, the heights of these peaks are suppressed in all 
three directions. 
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FIG. 4: Magnetization M versus the applied field H with 
D = (0.11, 0, 0) J (dashed line) and the pure Heisenberg model 
(solid line) for H//b. The corresponding shift Sq/2TT (o) is 
shown and also compared with the experimental data (-I-) 0- 

Figure Q shows the field dependence of the magne- 
tization M{H) and the momentum shift 5q for H//6. 
For comparison, the magnetization curve for the pure 
Heisenberg model and the experimental data for 6q are 
also shown in the figure. As seen from the figure, our 
numerical values of 5q agree quantitatively with the ex- 
perimental data 01 ■ However, we find that Eq. Q is 
not valid for copper benzoate. In particular, in contrast 
to M{H), Sq varies non-monotonically with H. We be- 
lieve this non-monotonic variation of Sq with H can be 
detected by neutron scattering measurements. In addi- 
tion, the magnetization is gradually saturated for large 
fields as a direct consequence of the DM interaction. This 
is different from in the pure Heisenberg model where the 
spins are fully polarized beyond a critical field (J-bHc = J, 
at which the two magnetization curves differ mostly. 
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FIG. 5: The staggered magnetization as a function of H for 
H//& and D = (0.11, 0,0) J 

When H//6, spins are polarized ferromagnetically 
along the 6— axis and antiferromagnetically along the 
other two directions. Figure ((SJ shows the staggered 
magnetization as a function of H for T) — (0.11, 0,0) J. 
At low fields, both M^„ and M^„ behave as H^/^. This 
field dependence of the staggered magnetization above 
2T deviates from the H^/'^\t}^^{J / ^bHs) behavior pre- 
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dieted by the efFeetive field theory |9| and is not sensitive 
to the D-vector. It may well be that the range of the va- 
lidity of the leading approximation in the effective field 
theory is narrower for the staggered magnetization than 
that for the gap. We should also mention that the numer- 
ical accuracy is lower when H is below 0.5T. In interme- 
diate fields, both M^,, and M^„ vary non-monotonically 
with H, in analogy with the field dependence of energy 
gaps. In the large H limit, AIs tends to approach zero. 
This means that the staggered field Hg induced by the 
external field is overwhelmed by the uniform field H in 
the high field limit, although Hg is approximately pro- 
portional to H. 
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FIG. 6: Local spin polarization for H = (0, 7, 0)r (left) and 
H = (0,25,0)r (right) for L = 120 and D = (0.11, 0,0) J. St 
is symmetric with respect to the middle point of chain, while 
the others are antisymmetric. Only is half of the chain shown. 

The above results show that the DM interaction affects 
strongly the properties of low-lying excitations although 
this term is much smaller than the Heisenberg exchange 
interaction. In particular, the staggered magnetic field 
induced by this term opens a gap in the spin excitations 
and leads to an antiferromagnetic long range order to 
accompany the ferromagnetic long range order induced 
by the applied field, shown in FigO The coexistence 
of orthogonal ferromagnetic and antiferromagnetic long 
range orders in an external field is a characteristic feature 
of quantum spin systems with DM interactions. 

In conclusion, we have shown that the microscopic 
model (1) including two inequivalent lattice sites and the 
DM interaction is a very good description of available 
experimental data. It can reproduce correctly the gap 
values in all three directions as well as the shift of the 
incommensurate peaks, without resorting to any addi- 
tional interaction like interchain coupling. The crossover 
to strong field behavior and the breakdown of the rela- 



tion between the momentum shift and the magnetization 
call for further experimental studies. Finally, we would 
point out that although the DM interaction is generally 
small, its effects could be very important for real materi- 
als. In particular, for metallic materials, e.g. Yb4As3 0|, 
one could use external fields to modulate the magnetic 
transport properties which is one of the focal points for 
the spintronics. 
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